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Abstract 

Classical oscillators of sextic and octic anharmonicities are solved analytically up to the linear power of A 
(Anharmonic Constant) by using Taylor series method. These solutions exhibit the presence of secular terms 
which are summed up for all orders. The frequency shifts of the oscillators for small anharmonic constants are 
obtained. It is found that the calculated shifts agree nicely with the available results to-date. The solutions for 
classical anharmonic oscillators are used to obtain the solutions corresponding to quantum anharmonic oscilla- 
tors by imposing fundamental commutation relations between position and momentum operators. 
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1 Introduction: 



The basic understanding of Physics is based on some physical models. The simple harmonic oscillator (SHO) is 
perhaps the most useful one among them. A particle subject to a restoring force proportional to its displacement 
gives rise to the model of a SHO. The Hamiltonian corresponding to a classical SHO of unit mass and unit frequency 
is given by 



and the Hamilton's equations are 



The equation of motion of the SHO is 



P 



dH _ 

dx 



(2) 



x + x = (3) 
where, equation (||) is used. Now the solution of the equation (|^) is found to be 

xo(t) = x(0)cost + x(0)sint. (4) 

The parameters x(0) and x(0) are the initial position and momentum of the oscillator. Thus the position of the 
oscillator at a later time t is completely known in terms of the initial position and momentum. Nevertheless, the 
momentum is also known from the Hamilton's equations (||). Hence, the oscillator problem is completely solved. 

The equation of motion of a quantum SHO may simply be obtained by imposing x and x as noncommuting 
operators. To avoid confusion, we represent the operators as capital letters instead of lowercase letters. Thus the 
noncommuting position and momentum operators obey the following relation 

[X(t),X(t)} = i (5) 

where h = 1. The solution of a quantum SHO in Heisenberg picture is the operator equivalent of the equation 
Now for real physical problems, the anharmonicity and/or damping are to be incorporated in the model 
Hamiltonian and hence in the equation of motion. However, the present letter takes care the anharmonicity only. 
The general form of the Hamiltonian of an anharmonic oscillator (having unit mass and unit frequency) of even 
order of anharmonicity is given by 

H = V — + — + —x 2m (6) 
2 2 2m y 1 

where m (> 2) is an integer and A is the anharmonic constant. Depending upon the problem of physical interests, 

different type of anharmonic oscillators will appear. For m = 2, the corresponding oscillator is called quartic 

anharmonic oscillator. Similarly, m = 3 and m = 4 give rise to the sextic and octic anharmonic oscillators 

respectively. Of course, the anharmonic constants A are different for different type of anharmonic oscillators. The 

equation of motion corresponding to the Hamiltonian (||) is 

x + x + Xx 2 " 1 - 1 = 0. (7) 
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where the equation (g) is used. 

2 Quartic Oscillator: Brief Review. 

The problem of classical and quantum anharmonic oscillators have already been studied by several authors. To 
begin with, we recall the problem of a classical quartic oscillator (Duffing oscillator). The equation (|?|), for m = 2, 
does not have exact solution. However, large number of approximate methods are available for the purpose of 
getting analytical solutions to the Duffing oscillator problem. These include perturbation technique Jlj], variation of 
parameters |§] and Taylor series approaches The ordinary perturbation technique, a pedestrian approach, leads 
to the unwanted secular terms 0. The removal of secular terms from the solution is a serious problem. There are 
some methods which are successfully used to summed up the secular terms for all orders. These include "tucking 
in technique" multiscale perturbation theory and renormalization technique jl). 

The search for the solutions of anharmonic oscillators enriched different brunches of Physics and Mathematics, 
for example, it is successfully used to enrich the subject of large order perturbation theory [^[jlo|], divergent expansion 
of quantum mechanics |^l| - |r^| , Laplace transformation representation of energy eigenvalues |fl4| and computational 
physics Hjl|. 

The solutions of quantum quartic anharmonic oscillators may be obtained in two different pictures. The 
Schroedinger picture deals with the time development of the wave function. On the other hand, the time de- 
velopment of operators are obtained in the Heisenberg picture. Basically, the Schroedinger picture is used to solve 
the quantum oscillators as eigenvalue problems. In these problems, energy eigenvalues are expressed as the sum 
of different orders of anharmonic constant A. The energy eigenvalues are found to diverge for large anharmonic 
constant. In case of small A, however, the eigenvalues for different orders are summed up and the convergence of 
these sums are ensured by the Borel summability [ |l3"| , [L9[ and/ or Stieljet conditions R], 

The c-number approach to the quantum oscillator problem is found very successful but the situation is different 
with the operator approach. Actually, the noncommuting nature of the operators (|J) pose a serious difficulty for 
the purpose of getting approximate solutions to the quantum anharmonic oscillator problems. Only few methods 
are available to-date. Aks pofland Aks and Caharat first obtained an operator solution of a quantum quartic 
anharmonic oscillator by using the method of Bogoliubov and Krylov. But the recent interest in the problem started 



with the work of Bender and Bettencourt p2f]23| who obtained the solution of a quantum quartic anharmonic os- 
cillator by using multiscale perturbation theory (MSPT). Later on Mandalj24|] proposed a Taylor series approach, 
Egusquiza and Basagoiti p"5|| developed renormalization technique, Kahn and Zarmi[p6| used a near identity trans- 
form method to solve the quantum quartic anharmonic oscillator problem. The Taylor series approach is found to 
be more successful compared to the other approaches for the purpose of getting second or higher order solution of 
a quantum quartic anharmonic oscillator p~7|| . 

In addition to the quartic oscillator, people have studied the higher anharmonic oscillators also JlH , p8|]30|| . 
These studies are made to obtain the eigenvalues as function of anharmonicity. More recently Pathak used a 
normal ordering technique to obtain a zeroth order multiscale perturbation theoretic solution oim — th anharmonic 
oscillator and Fernandez p2]| developed an eigen value approach to construct first order correction to the frequency 
operator for m — th anharmonic oscillator. However, a first order solution of higher anharmonic oscillators in the 
realm of Heisenberg approach are yet to come. Keeping that in mind, the present letter will take care the analytical 
solutions of classical and quantum oscillators of sextic and octic anharmonicities. These oscillators have potential 
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applications in the studies of nonlinear mechanics, molecular physics, quantum optics and in the field theory. This is 
the first time we give analytical solution of sextic and octic oscillators under the operator formalism. The solutions 
are used to obtain the frequency shifts of sextic and octic oscillators. The computed shifts are compared and 
found to have exact coincidence with the frequency shifts calculated by using a first order perturbation theory and 
frequency operators obtained by Pathak 

3 Sextic Oscillator. 

For m — 3, the Hamiltonian ^) and the equation of motion ([?]) correspond to the case of a sextic anharmonic 
oscillator. The solution for such an oscillator is obtained as the sum of different orders of anharmonicities. The 
corresponding solution is 

x(t) =x (t) +XX!(t) + ... (8) 

where, Xo(t) and X\{t) are zeroth and first order solutions respectively. The zeroth order solution (|]) is simply 
obtained for A = in equation (|7|). The purpose of the present section is to find the first order solution x\(t). The 
lower case x(t) and x(t) are to be treated as time developments of position and momentum of the classical sextic 
anharmonic oscillator. The equivalent operator representations are to be made by using upper case letters X(t) 
and X{t). Thus the solution of the classical sextic anharmonic oscillator is given by the following Taylor series 

t 2 

x{t) = x(0) + tx(0) + —2(0) + . . . (9) 

Here we assume that t is sufficiently small and a series (^) expansion is possible. Now we express higher order time 
derivatives of x{t) at t = as 

-2(0) - 5A.T 4 (0)i(0) 
x(0) + 6Ax 5 (0) - 20Ax 3 (0)± 2 (0) 

2(0) + 70Ax 4 (0)±(0) - 60Ai 3 (0)a; 2 (0) (10) 



x m (0) = 
x iv (0) = 
x v {0) 



where, we have neglected A 2 and higher order terms. We substitute the relations ( |l0| ) in equation (||) to collect the 
coefficients of x 5 (0), £ 4 (0)2(0), x 3 (0)x 2 (0), x 2 (0)i 3 (0), a;(0)i 4 (0) and i 5 (0) (without the constant multiplication 
factor A). Corresponding coefficients are 

C = -(1 - 6 + 71 - 1276 + 27741 - 656546+...) 

Ci = -(5 - 70+ 1275- 27740 + 656545-...) 

C 2 = -(20- 460+ 10680- 258920 + 6399340-...) 

(7 3 = - (60 - 1860 + 49320 - 12577720 + 31664580 + 

C 4 = -(120- 4320+ 120240- 3116640+ 78912360- 
and 

C 5 = -(120- 4320+ 120240- 3116640+ 78912360-...) 



(11) 



-...) 
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where the coefficient of x 5 l (0)i I (0) is denoted by Cj. The contributions directly from Taylor expansion are not 
taken into consideration. The r — th term of Co is given by 

tl = (-!)'' x 3^(25'' + 15x9'' + 240r- 16). 



Now we can write the r — th term of the coefficient of x 5 (0) in xi(t) as 



* r = (13) 
where the factor comes from the Taylor expansion part. The net coefficient of a; 5 (0) in x\(t) is obtained as 

K = ^2t r = ^( cos5i + 15cos3£- 16 cost - 120 tsint). (14) 

The parameter K exhibits fifth and third harmonic generations along with the secular term proportional to t sint. 
The secular terms are responsible for the divergent nature of the solution. This is an uncomfortable situation and 
has to be taken care properly in order to get a well behaved solution for all orders. Now, the remaining coefficients 
are obtained by using the same procedure as adopted for the evaluation of K . The corresponding coefficients are 

Ki = (5sin5t + 45sm3i - 280sint + UOtcost) 

K 2 = 3 | 3 (-5cos5f - WcosSt + 20cost - I20tsint) 

K 3 = ^(-bsinbt + 15sin3t - lAOsint + 120tcost) 

Ki = ^(5cos5t - AbcosZt + AQcost - \20tsint) 
and 

K$ = -^(sinbt — 15sin3t — 80sint + 120tcost) 

Hence, the first order solution is given by 



(15) 



x 1 (t) = ^2K i x 5 - i (0)x i (0). (16) 

i=0 

Thus the total solution (||) is simply 

5 

x(t) = x(0)cost + x(0) sint + X^KiX 5 ^ (0)^(0). (17) 



i=0 



Now to have a check of the solution (|l6|), we consider a special case x(0) — and x(0) = a. Hence the total 
solution reduces to an extraordinarily simple form as 

Aa 5 

x(t) = acost H (cos5t + 15 cosit — 16 cost — 120 tsint). (18) 

384 

It is clear that the presence of secular term proportional to tsint makes the life difficult as t grows. We confine 
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ourselves in the weak coupling regime and we assume C is a constant such that 



sin (CXt) = CXt 

■ (19) 

cos {CXt) = 1 



Now, the equation (|l8|) is rearranged as 



15 1 

x(t) = acos\t(l + — Aa 4 )l + Xa b ( cost H cosM H cosbt) (20) 

16 n v 24 128 384 ; v ; 



where the equation ( |19| ) is used. Thus the secular term present in the first order solution ( |18| ) is summed up for 
all orders and the frequency shift of the oscillator is obtained. The corresponding frequency shift of the sextic 
anharmonic oscillator is j^Xa 4 . The shifted frequency of the oscillator is also viewed as the renormalization of 
frequency by the anharmonic interaction. Interestingly, the frequency shift obtained in the solution (pp|) exactly 



coincides with that obtained by Dutt and Lakshmanan |3C||. Now to have an idea of how good the solutions (18) 
and ( po| ) are we compare them with the exact numerical solution obtained by Mathematica (figure 1). From figure 
1, we observe that the first order frequency renormalized solution ( po| ) coincides exactly with the exact numerical 
solution obtained by using Mathematica (for a = 2 and A = 0.01). However, the solution ( |l8| ) diverges with the 
increase of time (figure 1). The divergent nature of the solution manifests the presence of secular term. 




5 10 15 20 

t 

figure 1 

Now we replace x(t) and x(t) with the noncommuting operators X(t) and X(t) respectively. Finally, it is 
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straightforward to generalize the classical solution (||) into a quantum solution of the anharmonic oscillator by the 
imposition of commutation relation (H) . The corresponding solution of the quantum sextic oscillator is given by 



X(t) = X(0)cost + X(Q)sint + ^X 5 (0) x (cos5t + 15 cosM - 16 cost - 120 tsint) 

+ t4j[^ 4 (0)^(0) + X 3 (0)X(0)X(0) + X 2 (0)X{0)X 2 (0) + X{0)X(0)X 3 {0) + X(0)X 4 (0)} 

x (hsinht + 45sin3t - 280sint + \2Qtcost) 

+ iko I^ 3 (0)^ 2 (0) + X 2 (0)X 2 (0)X(Q) + X 2 (0)X(0)X(0)X(0) 

+ X(0)X 2 (0)X 2 (0) + X(0)X(0)X(p)X(0)X(0) + X(0)X(0)X 2 (Q)X(0) 

+ X 2 (0)X 3 {0) + X(0)X(0)X(0)X 2 {0) + X(0)X 2 (0)X(0)X(0) + X(0)X 3 (0)X(0)} 

x (-5cos5t - Ihcosit + 20cost - 120tsint) 

+ iku^ 3 (°) x2 (°) + X 2 (0)X 2 (0)X(0) + X 2 (0)X(Q)X(Q)X(0) + X(Q)X 2 (Q)X 2 (0) 

+ X(0)X(0)X(0)X(0)X(0) + X(0)X(0)X 2 (0)X(0) + X 2 (0)X 3 (0) 

+ X(0)X(0)X(0)X 2 (0) + X(0)X 2 (0)X(0)X(0) + X(0)X 3 (0)X(0)] 

x (— hsinht + 15sin3t - UOsint + 120tcost) 

+ 1535 [^(0)^ 4 (0) + X(0)X(0)X 3 (0) + X 2 (0)X(0)X 2 (0) + X 3 {0)X(0)X{0) + X 4 {Q)X(0)} 



(21) 



x (5cos5t - 45cos3t + AOcost - 120tsint) 

+ ^84° (sinbt — I5sin3t — 80sint + 120tcost) 



where X(0) and X(0) are the initial position and momentum operators. The solution ( |2l| ) may be written in the 
symmetrical form as 

X(t) = X(Q)cost + X(0)sint 

+ l±X 5 (0){cos5t + 15 cos3t~ 16 cost - 120 tsint) 

+ ^[X 4 (0)X(0)+1(0)X 4 (0)] 

x (5sin5t + Ahsinit - 280sint + \2Qtcost) 

+ ^[X 3 (0)X 2 (0)+X 2 {0)X 3 {Q)+3X(0)} 

x (— hcosht — Ihcosit + 20cost — 120tsint) 

+ ^ E [x 3 (o)x 2 (o) + x 2 {o)x 3 {o) + 3X(o)} (22> 

x (~5sin5t + lbsinM - HOsint + I20tcost) 

+ ^[x(o)i: 4 (o) + i: 4 (o)x(o)] 

x (5cos5t - 45cos3t + AQcost - 120tsint) 

+ 2A ^.g (0) [sinbt - 15sinM - 80sint + UOtcost) 

The equation ( p2|) is our desired solution for a quantum sextic anharmonic oscillator. Now we rearrange our derived 
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solution ( |22] ) with the help of (|l^ ) as 

X(t) = l{X(0)cos(t+^{H 2 + \})+co S (t+^{H 2 + \})X(0) 

+ X(0)sin (t + ^ {H 2 + i}) + , S zn (t + {H 2 + I}) X(0) 

+ 2A 384 °' > {cosht + 15 cos3t — 16 co.s£) 

+ ^[X 4 {0)X{0) +i:(0)X 4 (0)](5sm5t + 45sm3i- 280smi) 

+ ^[X 3 (0)X 2 (0)+X 2 (0)X 3 (0) + 3X(0)](-5cos5i-15cos3t + 20cost) (23) 

+ ^[1 3 (0)X 2 (0) + X 2 (0)Jf 3 (0) + 3X(0)](-5sm5* + Wsirtit- UOsint) 

+ 3^[X(0)X 4 (0) + X 4 (0)X(0)](5cos5t - 45cos3t + AOcost) 

+ 2 ^^4 (sinbt — Xbsinit — 80sint) \ 



where H = — ^ H j- 1 is the unperturbed Hamiltonian. Now the first order frequency operator, fi = 1 + 

X i H o + l) coincides exactly with the MSPT results of Pathak Q. The matrix element < n - l\X(t)\n > and 
hence the frequency shift of the oscillator may also be calculated with the help of following relations 

H \n> = (n+|)|n> 

X(0)|n> = ^ [Vn|n - 1 > +^/n + T\n + 1 >] forn ^ (24) 
X(0)|n> = ^ [Vn|n- 1 > -V^TT|n + 1 >] /orn 7^ 0. 

From the field theoretic point of view, n may be regarded as the number of photons present in the quantized 
radiation field. Now, the dipole moment matrix element in terms of the number eigenket \n > is 



<n-l\X(t)\n> = cos( ^ n) {< n - l\X{0)\n > cos (t + ffi {n 2 + 1}) 

+ cos (t+ ^ {n 2 + i}) <n-l|X(0)|n> 

+ < n- l|X(0)|n > sin(t+ Sf. {n 2 + |}) 

+ «n(t+5^{n 2 + |}) <n-l|X(0)|n> 



+ < n - 1 



2AX J (0) 
384 



(cos5i + 15 cos3t — 16 cos£) 



+ 3|i[^ 4 (0)X(0) + i:(0)X 4 (0)](5sm5i + 45sm3i- 280smi) (25) 

+ ^[X 3 {0)X 2 (0) + X 2 (0)X 3 (0) + 3X(0)](-5cos5t-15cosM+20cost) 

+ §^[X 3 (0)X 2 (0) + X 2 (0)X 3 (0) + 3X(0)](-5sin5t + lbsinM - UOsint) 

+ ^[^(O)^ 4 ^) + X 4 (0)X{0)](5cosM- AbcosZt + AOcost) 

, 2AX 5 (0) , 



'-(sin5t — I5sin3t — SOsint) \n > j 



Hence, we obtain the frequency shift of the oscillator as ^-(n 2 + 5). It is interesting to note that the frequency shift 
of the oscillator is possible even when there are no photons (i.e vacuum field) present in the radiation field. This 
phenomena is a direct outcome of pure quantum electrodynamic effect and has no classical analog. The frequency 
shift of the oscillator for a vacuum field is actually a second order effect and is reported elsewhere |27||. The solution 



(|23|) is of its first kind and does not have any direct check. However, the frequency shift in (|25|) may be compared 
with the frequency shift calculated by using first order perturbation technique [M. To have a check, we calculate 



S 



the energy of a sextic oscillator in n — th energy state which is given by 



1 5A 

En = + — (4n 3 + 6n 2 + 8n + 3). (26) 

2 48 

The energy difference between two consecutive levels is 

K\ 1 

A£ = £„-£„_! = 1 + — (n 2 + -). (27) 
Thus our calculated frequency shift and hence the solution agrees nicely. Finally we normalize the equation (|23f)with 



the help of the equation fl25[ ) as 



X(t) = 2cos ^ n) {X(0)cos (t + ^ {H§ + L}) +cos (t+™ {Hi + 1}) X(0) 

+ X(0)sin (t + ^ {Hi + i}) + sin (t + ^ {ff 2 + ±}) X(0) 

4- 2A 384 °" > {cosht + 15 cos3t — 16 cos/:) 

+ ^[X 4 (0)X(0) + X(0)X 4 (0)](5sm5i + 45sm3i- 280sm*) 

+ ^[X 3 {0)X 2 (0)+X 2 (0)X 3 (0) + 3X(0)]{-5cos5t~l5cos3t + 20cost) (28) 

+ ^[1 3 (0)X 2 (0) + X 2 (0)X 3 (0) + 3X{0)](-5sin5t + lbsin3t - UOsint) 

+ ^j[X(0)X A (0) + i: 4 (0)X(0)](5cos5t - 45cos3t + 40cosi) 



+ 



2 ^384^ (sin5t — 15sm3i — 80smi)| 



The equation (|28| ) is the desired solution for quantum sextic anharmonic oscillator where the secular terms for all 
orders are summed up. 

4 Octic Oscillator. 

Depending upon the nature of nonlinearity, the higher anharmonic oscillators will come into picture. For example, 
octic anharmonic oscillator may be obtained if we put m — 4 in the equations (^|) and ([?]) respectively. Using the 
similar procedure (as adopted for sextic oscillator) , the solution (up to the linear power of A) for an octic oscillator 
follows immediately as 

x{t) = x(0)cost + x(0)sint 

+ (coslt + Ucos5t + 126cos3i - Ulcost - 8A0tsint) 

+ AA( 3072 (0) (7sin7t + 70sin5t + 378sint - 2373sint + MQtcost) 



3072 

X±2 f^ {0) {-21cos7t - l26cos5t - 126cos3i + 273cosi - 2520tsint) 
A ^3072 4(0) (-35sm7i - 70sin5t + 630sin3t - 3815sint + 2520tcost) 
A ^ 4 3072 3(0) (35cog7^ - 70cos5t - 630cos3t + mhcost - 2520tsint) 
X±5 f^ (0) (21sin7t - 126sin5t + I26sin3t - 2415sint + 2520tcost) 
^3072 (0) (~7 c os7t + 70cos5t - 378cos3t + 315cost - MOtsint) 
¥^ {-sin7t + lAsinbt - 126sm3i - 525sint + MOtcost) 



(29) 
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To check the validity of our solution, we consider a special case ar(0) = b and i(0) = 0. After a little algebra we 
have 

x(t) = bcos(l + ^b 6 )t + ^—(-Ulcost + 126cos3t + Ucosht + coslt). (30) 
128 3072 

Thus the frequency shift of the classical octic oscillator is proportional to the sixth power of amplitude as long as 



the first order solution is concerned. The solution (30) as well as the frequency shift have exact coincidence with the 



solution obtained by the procedure introduced by Bradbury and Brintzenhoff f33|. Now we generalize the solution 



(29) to have the corresponding quantum solution in the symmetrical form as 



X(t) = X(0)cost + X(Q)sint 

+ x *^° ] (coslt + Ucosbt + \26cos3t - Ulcost - MOtsint) 

+ sm[X(0)X 6 (0) + X 6 (0)X(0)](7sin7t + 70sin5t + 378sint - 2373sint + 840tcost) 

+ ^[X 5 (Q)X 2 {0) + Z 2 (0)X 5 (0) + 10A 3 (0)} 

x (~21cos7t - \26cos5t - \26cos3t + 27'icost - 2520tsint) 

+ gA_[X 4 (0)X 3 (0) + A 3 (0)A 4 (0) + 9(A 2 (0)A(0) + A(0)A 2 (0))] 

x (-35sin7t - 70sin5t + 630sin3t - 3815sint + 2520tcost) (31) 

+ gA_[X 3 (0)A 4 (0) + A 4 (0)A 3 (0) + 9(A V2 (0)A(0) + X{Q)X 2 (0))} 

x (35cos7t - 70cos5t - 630cos3t + 665cost - 2520tsint) 

+ ^[X 5 (0)X 2 (0) + X 2 (0)X 5 (0) + 10A 3 (0)] 

x (21sin7t - 126sin5t + 126sin3t - 2415sini + 2520tcost) 

+ 3(^2 [X 6 (0)X(0) + X(0)X 6 (0)](-7cos7t + 70cos5t - 378cos3t + 315cost - MOtsint) 

+ A f 7 ( 2 0) (-sin7t + UsinM - 126sin3t - 525sint + MOtcost) 

The problem of secular terms are taken care by the same procedure as it is done in the case of a sextic anharmonic 
oscillator. The final solution appears as 

X(t) = - — r 35xL 2 , a wi (x(0)cosD,t + cosSltX(0) + X(0)sinSlt + sinCltX(0) 

2cos[ (on -\-o)t\ ^ 

+ 2A 3 ^ 0) (cos7t + Ucosht + 126cos3t - Ulcost) 

+ [X(0)X 6 (0) + X 6 (G)X(0)] (7sin7t + 70sin5t + 378sint - 2373sint) 

+ [A 5 (0)A 2 (0) + X 2 (0)X 5 (0) + 10A 3 (0)] 

x (-21cos7f - 126cos5t - \2Qcos3t + 273cost) 

+ Mf2i x4 (0)X 3 (0) + A 3 (0)A 4 (0) + 9(A 2 (0)A(0) + X(0)X 2 {0))} 

x (— 35sin7t — 70sin5t + 630sin3t — 3815sint) 

+ 3 ? ^[A 3 (0)A 4 (0) + A 4 (0)A V3 (0) + 9(A >2 (0)A(0) + A(0)A 2 (0))] 

x (35cos7t — 70cos5t — 630cos3t + 665cost) 

+ [A 5 (0)A 2 (0) + X 2 (Q)X 5 {0) + 1QX 3 (0)} 

x (2\sin7t - \2Qsinbt + 126sin3t - 2415smt) 

+ Mr! [^ 6 (0)A(0) + X(0)A 6 (0)](-7cos7i + 70cos5t - 378cos3t + 3lbcost) 
xx 

3072 



2XA < sin7t + Usin5t - I26sin3t - 525smt)| 
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where Q = 1 + (4i?g + 5-ffo) , is the first order frequency operator which is in exact accordance with the frequency 
operator obtained by Pathak[ 

]. Now, the matrix element < n— l\X(t)\n > is calculated to obtain the shifted frequency as a/ = l + ^(n 3 +2n). 
Interestingly, the frequency shift 3 j£(n 3 + 2n) depends on n 3 . It is clear that the vacuum field has no role to play in 
the frequency shift of the quantum octic oscillator. This is the first time we give the solution of an octic oscillator 



( |32[ ) and therefore, it is not possible to have a direct verification of the solution. However, to have a possible 
verification of the solution, we calculate the frequency shift by using Rayleigh-Schroedinger perturbation theory. A 
first order calculation under this theory shows that the energy of the octic oscillator in n — th state is given by 

E n = (n + ~) + ^-(^ + An + 5n 2 + 2n 3 + n A ). (33) 
2 o4 2 

The energy difference between two consecutive levels is AE — E n — E n -\ = 1 + ^-{n 3 + 2n). Thus, the frequency 
shift exactly coincides with the frequency shift obtained by us. 



5 Remarks and Conclusions: 

We obtain analytical solutions (first order) for sextic (||) and octic j2g| ) oscillators in classical pictures by using 
Taylor series method. The solutions exhibit the presence of usual secular terms. The secular terms for all orders are 
tucked in to obtain the frequency shifts of those oscillators. It is found that the calculated shifts agree satisfactorily 
with the shifts obtained by other methods. 

The position and momentum of classical oscillators are replaced by their corresponding operators and the 
commutation relation is imposed. Hence, the solutions for quantum sextic and octic oscillators in the realm of 
Heisenberg picture is obtained. Interestingly, the vacuum field has no role to play in the frequency shift of the 
quantum oscillators of quartic and octic anharmonicities. However, the vacuum field causes the frequency shift of 
a quantum sextic oscillator. These apparent discrepancies are easily understood in terms of field theoretic point of 
view where the vacuum field effect is considered as a second order one. In fact, the frequency shift due to vacuum 
field is obtained in a second order solution (i.e solution contains A 2 ) of a quantum quartic oscillator p7[ . 

To the best of our knowledge, there only few possible methods available for the purpose of getting analytical 
solutions of a quantum quartic anharmonic oscillators |20[ - p7| |. The method devised by one of us [Q is already 
tested its usefulness for the purpose of getting an analytical second order solution of a quantum quartic anharmonic 
oscillator j27j]. However, there is no publication till date which deal with the solution of higher anharmonic oscillators. 
Therefore, the solutions for higher anharmonic oscillators of the present letter will serve a lot of academic interest. 
In addition to these the solutions will be useful in nonlinear mechanics, molecular physics and in the field theory. 
In particular, the solutions for quantum anharmonic oscillators are useful to investigate the quantum statistical 
properties of radiation field. For example, the squeezing effects are reported for a third order nonlinear non- 
absorbing medium by using the solution of the quantum quartic oscillator ||27| , |35|l . 

The present solutions along with the solutions given by us seem to be useful for the purpose of getting 

a general solution of m — th anharmonicities. Such solutions will be taken care in future publications. 
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